Modification of Eckart theory of relativistic dissipative fluid 
by introducing extended matching conditions 
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We deal with a novel approach to formulation of the relativistic dissipative hydrodynamics by 
extending the so-called matching conditions widely used in the literature. The form of the non- 
equilibrium entropy current can be determined by requiring thermodynamical stability of the entropy 
current under extended matching conditions. We derive equations of motion for the relativistic 
dissipative fluid based on the Eckart theory and show that linearized equations obtained from them 
are stable against small perturbations. It is also shown that the required fluid stability conditions 
are related to the causality of the model. 
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I. INTRODUCTION 

Israel and Stewart (IS) provide phenomenological 
framework for a relativistic dissipative fluid [11, 
In their model, a possible general form for the 
non-equilibrium entropy current is described by the 
dissipative part of the energy-momentum tensor and by 
the particle current up to the second order in the devia- 
tion from the equilibrium state. After the stability and 
causality of the IS theory have been shown by Hiscock 
and Lindblom the causal dissipative hydrodynam- 

ical model was adapted to study the dynamics of hot 
matter produced in ultra-relativistic heavy-ion collisions 
by Muronga [J 01- In recent years, the IS model plays 
important role in the analysis of the experimental data 
obtained by RHIC and LHC (see for example [1]). 

In parallel to application of IS theory, investigations of 
the basis of relativistic dissipative hydrodynamic theory 
was vigorously continued The reason is that the 

IS theory in its present form is too general and complex 
from the point of view of the quantum chromo dynamics 
(QCD), which is belived to be the fundamental theory 
for strongly interacting system Q. Another reason is 
that theory of the relativistic dissipative hydrodynamics 
is not yet fully understood because, for example, the 
equation of motion of the fluid it uses depends on the 
choice of the Lorentz frame (or on the definition of 
the hydrodynamical flow [l^). Since the dissipative 
part of the energy-momentum tensor ^T'"^ and the 
particle current SN^^ cannot be determined uniquely 
by the second law of thermodynamics, one usually 
introduces some constraints to fix them, they are known 
as matching conditions. 

The matching (or fitting) conditions are also intro- 
duced because of necessity of matching the energy density 



and baryonic charge density (e, n) in a non-equilibrium 
state to the corresponding equilibrium densities: (Seq, 



j), e = £oq, n — "•oqj or equivalently, 



Su^N^' = 0. 



(1) 



These matching conditions allows to determine the ther- 
modynamical pressure Peq(eeq, f^eq) (defined as work 
done in isentropic expansion) via equation of state for the 
equilibrium state. Here the Pcq should be distinguished 
from the bulk viscous contribution, 11 = — ^Afj_i, ST '^'^ , 
present in the energy- momentum tensor Q • The match- 
ing conditions are also needed because they are necessary 
for the thermodynamical stability of the entropy current 
(see Appendix A in Ref.flSj). However, the matching 
condition given by eg. (IT]) are not unique. So far, except 
of some recent works |l4l . [l5j , they were not investigated 
in detail. 

A state of relativistic dissipative fluid is described 
by energy- momentum tensor T'^'^{x) and by the baryon 
number current (x) , which obey the conservation laws 

M 



= 0, 
0, 



and the second law of the thermodynamics 



(2a) 
(2b) 

(2c) 



Because of the uncertainty in definition of the fiow veloc- 
ity u^{x) for a non-equilibrium fluid, one needs, unlike in 
the case of perfect fluid, to flx the frame for the fluid con- 
sidered. Throughout this paper we employ the so-called 
Eckart frame |20| . It means that the hydrodynamical flow 
velocity (with normalization u^Ufi = 1) is defined by 
using baryon charge current 



u^(a;) 



ATM 



(3) 
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In this frame one has always A'tN^ 



with A'' 



5^ — u^Ui, being the projection operator orthogonal to 



2 



the four vector u^^. Ahhough the alternative Landau- 



frame 21| (in which the nonzero baryon number cur- 
rent exists A^iV'*' 7^ but the energy flow disappears, 
= Uy/S.'^T^^ — 0) is more relevant in the context of 
ultra-relativistic heavy-ion collisions, we use the Eckart 
frame because of our previous work 22| , in which we have 
argued that for the relativistic fluid in the presence of the 
long-range correlations it^w^^T'"^ ^ 0. 

In the present work we reformulate the relativistic dis- 
sipative hydrodynamics based on Eckart prescription j20| 
using following extended matching conditions, 



i^u^STf^" = 0(1), u^Sm ^ 0(2), 



(4) 



instead of those given by eq.([T]). Here, 0(1) denotes the 
first order infinitesimal quantities, i.e., the difference of 
fiuid dynamical quantities between the equilibrium and 
the dissipative state (for example, bulk pressure 11 = P — 
Peq is of the 0(1) order). On the other hand, u^SN^^ is 
assumed to be of the order 0(2) because of introduction 
of infinitesimal 0(2) to the condition u^^SN^ — 0; it is 
sufficient to derive the equations of motion for stable fluid 
as discussed in IIIII 

Hence, the main purpose of this article is to investigate 
the system described by equations eq. ((2a|) and eq. (j2bl) . 
in the frame defined by eq.(I3]), while observing eg. (Pel) , 
with constraints provided by ^ instead of eq. ([1]) . 

This article is organized as follows. In Sec|lTl we 
briefly review the relativistic fluid dynamics with stan- 
dard matching conditions. Then extended matching con- 
ditions are introduced into the hydrodynamical model 
using general form of an irreversible entropy current. In 
Sec lIIIi we check the stability and causality of the fluid 
obtained from our model. We close with Sec. IIVI con- 
taining the summary and some further discussion. Some 
technical details are presented in Appendix A and B. 



where T^^ and iV^q are equilibrium energy-momentum 
tensor and baryonic charge current: 



eq 



(7a) 
(7b) 



In the equilibrium case, the energy-momentum conser- 
vations T^f^ = 0, and the baryon number conservation 
^cq-fj. — Oj together with thermodynamic relations, result 
in the locally conserved entropy current: 

^o^q;^ = -«,M^c'; + PX;,T,\^ + W'Poq] = 0. (8) 

To this current we introduce now dissipative corrections 
by adding corresponding dissipative corrections STf^" and 
SN^ to the energy-momentum tensor and to baryon num- 
ber current appearing in (jj^ and ^q'^ in eq.([6]). In this 
way one extends the expression of equilibrium entropy 
current towards the off-equilibrium entropy current. 



X^A„^<5r"^. (9b) 

Notice that, because of the thcrmodynamical stability, 
term proportional to A^^iJT"'^, which appears due to the 
natural extension of ea. (j6b[) . is usually dropped. How- 
ever, since we shall generalize these conditions, we re- 
tain this term and introduce a phenomenological param- 
eter X which determines its contribution to the non- 
equilibrium entropy current. Note that in the Eckart 
frame bN^" = 0. This is because A^(57V^ = (by def- 
inition of the Eckart frame) and because of one of the 
matching conditions, u^JiV^ = 0. From another match- 
ing conditions, u^Uiyi5T'^^=0, it follows that one can also 
drop term proportional to u^u^ . As result one gets 



-HA^"^ + M^^u" -I- Wu^" 



(10) 



II. RELATIVISTIC DISSIPATIVE 
HYDRODYNAMICS WITH EXTENDED 
MATCHING CONDITIONS 



Eckart theory with standard matching 
conditions 



The entropy current can be written in the following 
simple form 



(5) 



where 0^ and ^^'^ are, respectively, vector and 2-rank 
symmetric tensor, a = /i&/r and — I3u^ with /3 = 
1/T, T and /ifc being, respectively, the temperature and 
baryon chemical potential. In a local equilibrium case it 
is then given by 



^cq 



ATM 
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2 ^apT^q, 



(6a) 
(6b) 



where H is bulk pressure, is energy flow and tt^'' 
is shear tensor. Using now this expression for 5T^^ in 
eq. (|9bp we obtain following possible expression for the 
non-equilibrium entropy current 



H , 



(11) 



In the usual formulation, i.e., for standard matching con- 
ditions, term proportional to Hu^ is dropped to prevent 
the occurrence of instability in the entropy current (in 
our case it would mean assuming x = 0). In this case 
one gets the known first order theory form of the entropy 
current 12011: 



(12) 



Correspondingly, the entropy production in the Eckart's 
1st order theory is given by 



3 



However, for more general form of matching conditions 
proposed here this term should be included with param- 
eter X used to restrict the form of the entropy current in 
the non-equilibrium state. 



B. Extended matching conditions 

In what follows, we propose and discuss in detail the 
extended matching conditions, 



0{2) infinitesimal quantities. In the most simple way 
they can be then written as: 



A = K.n, 

/ib(5n — - 



(20a) 
(20b) 



u^uxST^^' - A, 
UaSN" = 6n, 



(14a) 
(14b) 



to be used instead of the standard matching conditions 
given by eq.([T]) and widely used in the literature. In this 
case the dissipative part of energy-momentum tensor is 
give by 



where k,^,^', and ^" are positive constants. Note that 
and ^" have dimension [GeV]^^ whereas k has is 
dimensionless. It is interesting that A — 0(1) (k = 3) 
is obtained in Ref. [ISJ and also obtained independently 
in Ref. m (see also Ref. (Hi). 5n = 0(2) is needed 
to stabilize the so-called transverse mode propagation of 
the dissipative fluid (discussed below in Scc lIII A]) . 

A comment concerning the matching conditions may 
be in order here. The IS theory introduces contributions 
of the second order, 0(2), Q'^((5r^V^^) to the entropy 
current (see eq.(2.17) in the ref.i2]). In the present 
model, the second order contributions to the entropy 
current are brought in by the matching condition 



5T^"' = Ku^u" - + W^u'' + W'u^ + 71^" (15) eg. ((20b)) . u^5N^' ^ (together with contributions of 



instead of eq. (jlOp . Correspondingly, the entropy current 
eq. ([5a|) is modified and takes form 



't 



S^^ + — - aSnu^" + Kp^" - xl^^u^" . (16) 



If we now require that 

-aSnu^" + A^f" - x/SHm^ = 0, 

or equivalently, that 

-^b5n + A-xn = 0, (17) 

the entropy current reduces to its usual Eckart's form, 
I— S'^. Eq. lfTTl) defines therefore new thermodynani- 
ical stability condition which replaces previous require- 
ment that in eq. ljlip X — 0; this new condition does not 
affect the local entropy density. On the other hand, ex- 
tended matching condition, eq. (|14ap . results in 

therefore now one gets 

^^ = 5c^q + -^ ^S^^ + 13,ST^''-A13'^. (18) 
The entropy production is thus in this case given by 



the first order, brought by u^u^^T^" ^ 0, cg. pla)) . 
Notice, however, that, because of the stability condition 
eq. p?)) . one has = 0, i.e., terms introduced by 
extended matching conditions do not contribute to the 
total entropy current. In this way our model is, in fact, 
a first order model and because of this differs from the 
IS approach. 

Notwithstanding the fact that out matching conditions 
do not give direct contribution to the entropy current, 
the resulting equation of motion of fiuid changes its form 
accordingly. The thermodynamical stability condition, 
eq. p?)) . now reads 



X 



en - e 



n 



e^^^. (21) 



Because the entropy production eq. (|T9|) can be rewritten 
as 



— (f^bSn)[——] - — 

dT fj,h dr dr 



Xm, (22) 



the second law of thermodynamics is guaranteed (with 
C, A and rj being, respectively, the bulk viscosity, heat 
conductivity and shear viscosity - all positive constants). 



2t] 



> 0, 



(23) 



dSn dA 



[aSn - liA]e, (19) 



where 9 = denoting the local expansion rate. The 
last two terms represent corrections to the usual Eckart's 
formula, eq. (IT3t . 

We assume now that A and 5n can be expressed by the 
respective scalar quantities: 11 and W^Wf^ and n^^" n ^j,^ . 
Furthermore, we assume that they are given by 0(1) and 



provide that the following constitutive equations are 
hold: 

7 =-(l + x)^^-(^ + 20f^ + (^f^^)n, (24a) 
Q 11 dr /ih dr 



XT 



-u'u,.^,-2e^ + i^^)W„ (24b) 



T 



^ dT 



dr ^/ib dr 
fib dr 



(24c) 
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Because eq. (|24ap includes term proportional to 
1/n dn/dT, one can introduce for the bulk pres- 
sure n an arbitrary constant z (with dimension [GeV]^) 
and write: 



[l-^n]-+2a- + In- 

z dr z z dr z 



c 



{1 + X)0, (25) 



where 7n^C^^ 



C. Small perturbations of the dissipative fields 

To check the stability and causality of our model, 
we shall consider small perturbations of H, and 
TT^'' fields. The bulk pressure 11 can be written as 
n = Ho + 5Ii with the background reference field XIq 
and its small perturbation field STL field [l^. Identi- 
fying arbitrary constant z with Hq and noticing that 

ln(l -I- IP) sa ^'^n, one can rewrite the above equa- 
tion as the following a set of equations: 



dsn 



r 1 + X 



+ [1 - 7n]OT = -C(l + x)5e, 
where the relaxation time rn is given by 
rn = C (2^ + «/no) , 



(26a) 
(26b) 

(27) 



and 9 = + 59. The ea. (l26ap represents a background 
Ho field and eq. (j26b[) represents equation of motion for 
the perturbed field of the bulk pressure 511 25 1. 

Similar approach can be also applied to W^^ and tt'"^, 
with perturbation of the heat flow and the shear viscosity 
around — and tt'^'^ ~ 0, leading to 



dSW^ 

dSnt^" 
dr 



- [1 - lw]SW^ = AT 
[1 - 7^](57r^'' = 277V 



. T 



A u 

-u < 



(28) 
(29) 



Here r^, and denote the respective relaxation times: 



/jb ar ' ' " fit, dr 

For simplicity, we drop in what follows terms with 7n , 
7u, and (assuming that ^^§7- ~ 0) and arrive at the 
following constitutive equations: 



dSWf' 

d5^T^"' 
dr 



= AT 



rV^T 

(as, I') 



A M 



(30a) 
(30b) 
(30c) 



III. STABILITY OF THE FLUID OBTAINED 
FROM OUR MODEL 

The stability of general class of dissipative relativistic 
fluid theories was investigated by Hiscock and Lindblom 
Denoting by SV{x) the difference between the ac- 
tual non-equilibrium value of a fleld V{x) and the value 
in the background reference state, Vo{x), we assume that 
variations SV are small enough so that their evolution is 
adequately described by the linearized equations of mo- 
tion describing the background state. We shall now in- 
vestigate the stability of the fluid obtained in our model 
following prescription proposed in Ref. In what 

follows: 

1. The background reference state is assumed to be 
homogeneous in space. Notice that, unlike in Ref. 

in our case it is not an equilibrium state but 
rather a non-equilibrium one with 11 = IIq and 
with = TT^" = 0. Furthermore, the background 
space-time is assumed to be flat Minkowski space, 
so that all background fleld variables have vanish- 
ing gradients. 

2. We consider following plane wave form of pertur- 
bation propagating in x direction 



5V = (5Vo exp(?;A:a; + Fr) 



(31) 



Linearized equations for dissipative fluid dynamical 
model are given by 



5[Tn■.^^ = 0, 



(32a) 
(32b) 



with the perturbed energy-momentum tensor and baryon 
number current given by, respectively: 

+ (fcjq + K6Tl)u^'u'' - {SP:^ + OT)A^^ 
+ SWu" + SWu'' + St:'"' , (33a) 
SIN''] = Sncqu'' + noq(5u^. (33b) 

Here £*q and P*^ are energy density and pressure in the 
background non-equilibrium state 

^cq = £oq + kHq , P*q = Pcq + Ho . 

However, since Hq = 0{1) and (5[Ho] — (it has vanish- 
ing gradient and is constant in r), terms proportional to 
Hg do not contribute to the linearized equations eq. (j33ap 
(with only terms up to the order 0(1) present). Hence, 
one can replace in the ea. p3ap e* and P* by the, respec- 
tively, Eoq and Pcq. 

The perturbed fluid dynamical flelds must satisfy con- 
straints 

u^Suf" = 0, SWfiu" = 0, Stt^^u^ = 0. 

Hence, in the rest frame of fluid, the proper time t 
component of the flow velocity fleld vanishes, Su'^ = 0. 



5 



Moreover, since the contribution from 5n is C(2), the 
term Snu'^ appearing in eq. P3bp can be disregarded. We 
obtain therefore foUowing hnearizing equations for the 
energy- momentum tensor and baryon number current: 

+ (r&eq + Krsn)u'' - (V^Poq + vm) 

+ {ikSW)u'' + rSW + {ik)STT^'' = 0, (34) 
S[N''lf,^TSn,q + n,q{ik6u'=) = 0. (35) 

The hnearized constitutive equation for 11 have now the 
form (with k = 1 + n ) : 



where the matrices Q and R are given by 

Q = 
/ 






r 


ikfl^q 








r 





ikhcq 


kT 


ik 


zk^ 


^k^ 


r/leq 


ik 


r 








ikk 


l+rnr 

c 







ik dT 


ik dT 


r 


i+r„r 


T Oeoq 


T Oricq 


AT 








ik 









R 



hcqV r 
r 

ik 



XT 





ik 






ik 



4^/3 



(42a) 
(42b) 



(1 + mV) 
C 



(5n = -R{ik)5u^, 



(36) 



whereas for and tt^'', we have 



(l + r„r)^^^^ V^<5T dSu^^ 



XT 

(1 + r^r) 



27] 



T dT 



(37) 



(38) 



The parameters ^, and ^" introduced in eq. (j20bp 
have been absorbed in the expressions for relaxation 
time, Tn,Tu, and Ttt, respectively. On the other hand 
parameter k, in the expression of the entropy production 
is kept not absorbed in rn. Its role will be discussed later. 

All perturbation equations can be expressed in concise 
matrix form: 



MpY^ = 0, 



(39) 



where 5Y^ represents the list of fields. The system ma- 
trix can be expressed in a block-diagonal form when 
one chooses the following set of perturbation variables [4] 



5Y^ = { fccq, (Jncq, 8u^, m, SW, (Jtt^^, 



(40) 



In this case. 



Q 



M 



R 



R 



(41) 



respectively, and I is the 2x2 unit matrix. The 
denotes the enthalpy density which is defined by 



"cq 
^cq — ^cq 



P. 



cq- 



For r and k satisfying dispersion relation 

[detM] = [det R]2[det Q] = (43) 



one has plane- wave solution such like eq.([3T|) for the lin- 
earized equations of the system eq. . 

The set of all exponential plane- wave solutions is char- 
acterized by the collection of roots obtained by setting 
the determinants of matrix R and Q separately equal to 
zero 

In what follows we shall discuss in detail the stability of 
transverse and longitudinal modes separately and provide 
the general stability conditions required in our case. 



A. Transverse mode 

The dispersion relation obtained by setting 

[7?][AT]dct(R) 



3 

?n— 







(44) 



corresponds to the solution of the perturbation equation 
which is referred to as the so-called transverse mode. The 
coefficients r™ are given by 

Tj, = Tw [hcqTw - XT] , 
r2 = hcq{Tu, +T„) - XT, 

ri = heq + r]T^k'^, 
ro = rjk"^. 



When all coefficients of the 3rd order equation (j44|) have 
the same sign, the three solutions consists of a pure real 
root with negative sign and two complex number roots 
with negative in real parts. In this case, since each real 
part of those solutions is negative, the general solution 
which is linear combination of those three solutions, is 
stable. Since rg and ri are positive defined, then the 



6 



stability condition sought after is that and r2 must be 
simuhaneously positive, i.e., that, 



AT 



— < /icq and 



< h. 



cq- 



(46a) 



Then, the necessary and sufficient condition which simul- 
taneously satisfies ea. (|46ap is 



— < h, 



cq ■ 



(46b) 



When condition eq.(j46b| is satisfied, the transverse mode 
solution of the dissipative fiuid dynamical equations can 
evolves without facing an instability problem. 

Note that for rn — Tu, — t,^ — (with also k, ^, ^' and 
are all zero but with C^^^V 0, i.e., in the case of 
the standard matching condition), one gets the disper- 
sion relation obtained already by Hiscock and Lindblom 

in Ref.g, T = ^ ± \J + It has two real 
solutions the linear combination of which leads to well 
known instability of the transverse mode. 



B. Longitudinal mode 

The determinant of the coefficient matrix Q can be 
expanded in the following way: 



det(Q) = (ikV) det(Qi) 



[i±^];.(^fc) det(Q2) 
.l + r^r„l + Tnr 



477/3 



where 



Qi = 



ik 



Q 



2 — 



Q3 



r 

c?eoq 



ik dT 
T aeoq 



r 

ik dT 
T aecq 



r 

%k~^ — ~ 

oeoq 
1 ik dT 
V T Seoq 



ik 

l-l-rnr 
C 



r 


9ncq 
ik dT 
T Oncq 

r 


ik dT_ 

cq 



c 



ik 

r 




] det(Q3), (47) 



T dn, 



XT 


nV 
ik 


ikric 
ikhc 
The, 

r 



ik 







ik 

r 

AT / 

q 
q *fc 

, r 

AT 



(48a) 



(48b) 



Then the frequencies of the so-called longitudinal mode 
are given by the roots of the following dispersion relation: 



4?7 AT C 
3r^ Tn 



]det(Q) 



6 

E 

n=0 



g„r" = 0, 



(49) 



where the coefficients g„ are given by 
AT . 



96 



95 



XT 2 



3h 



cq 



94 = 



93 = 



T 

AT 



— a 

m 

XT - 3fh. 



— [k—X - J2) + - — 



Jl 



cq 



kC 2 

Tn T^TT 
An 2 

'J'^TT 'Hilt; 



XT^J^CX 2J2^ 

Tu, Tn Ttt TnTT 



3Ji 
f 



^2 _ ^cq 



AT 



92 



91 = 



—X 

m 



XT J2 



XT 



kC X 



Tn Ttt 



4 77 

3 Tjr 

4?7 1 

STtt TnT„ 

4?7 X 
3t„ Tn 



3 ' 

3f 

2J3 
TnTT 



a; 



+ 



Jl 



Tn T^T^ 



9o 



AT J, 



Tu; TnT^r 



The notations used for fa.b(for a, 6 = 11, w and tt), f and 
T is: 



Tah 2 Tq 

r=-(Tn 



1 , 1 1,1 1 1 , 

n T 6 Tn Tyj Tjr 



(TnTtuT^ 



sl/3 



whereas J„ (ri = 1, 2, 3) are defined by 



J , 5Peq 
•Jl = ^cq-:: 

oe, 

J 5Peq 



ap 



oq 



''cq 



cq 



''cq 



dn. 
dT 



cq 



dec 



T dn^ 



J "-cq T^-^ cq 

J?, 



dPra dT 



dT dR 



T l9eeq arieq 



l9£eq l9neq 



The definitions of x and X are 



1 



dP, 



cq 



dSc 



X 



1 9T 

Tdz: 



Notice that t (and quantities such as f and t associated 
with it) are positive. Also X is positive, but the sign of 
X depends on and n. In what concerns the signs of 
Jl, J2 and J3 it occurs that, apparently, Ji > 0. On the 
other hand, because I < and other three factors 
in J3 are positive, then we have J3 < 0. In what concerns 
J2, since one has (see Appendix iBjl 



J2 = ^\s. 



dr 



cq 



dn, 



cq 



cq 



dfl 



^cq 



dT 



^cq ^S^cq 1 

~~djr^' 

^l^"^-,^'' > 0, we have J2 < 
0. To observe J2 < 0, consider J2 in the Boltzmann 



where J is Jacobian J 
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approximation. In the limit, we have 



J2 = ^ 



1 (^cq £^cq)^cq 



J y2 

which has negative sign apparently. 

C. General stability conditions 

To get stable evolution of the fluid, both transverse and 
longitudinal mode must be simultaneously stable against 
small perturbation and each condition needs to be satis- 
fied in a without a contradictory. The longitudinal mode 
needs to satisfy the requirement that, as the same as the 
transverse mode, all coefficients g„ in ea. (|49)) have the 
same sign in order to obtain solutions of plane- wave type 
ea. (l5T]) with negative real part of complex F, Re F < 0. 
Since go (and also qi) is negative then all the coefficients 
qi must be negative, < (z = 2, 3, • • • , 6). Now condi- 
tions that < 0, < and the second term of are 
negative, are equivalent to demanding that respectively. 



AT 



< K 



XT 



<[1 + 



1/Tu 



and 



— < IH \hc 



The necessary and sufficient condition for all these three 
quantities to be negative is that 



AT 



< K 



(51a) 



This is exactly the same condition as that for the sta- 
bility of the transverse mode, cf. eq.(j46b|. To continue, 
conditions that the first term of 54, the 93 and q2 to be 
all negative is that a; > 0, i.e., that 



cq 



de, 



cq 



1 dP, 



cq 



T ds 



cq 



< 



(51b) 



where = I is adiabatic velocity of sound. From 
eq. (|51bl) it is clear that parameter k, is related to the 
speed of sound wave, Cs, propagating in the fluid. Note 
that, when k — > in eq. one flnds that bound on 

the speed of sound, Cs < 1, can be exceeded violating 
causality. On the other hand Cs > 0. Therefore eq. (j51bp 
gives the following restriction on the parameter n: 



1 dP, 



cq 



1 

< - < 1 
K 



1 dP, 



cq 



T dScq K T dScq 



(52) 



To summarize this part, the stability condition for a 
dissipative fluid in an off-equilibrium state with pressure 
P = Pcq + Ho (the corresponding energy density is given 
by £ = Ecq + kYIq) is that 



< Ecq + Pc 



- cq 



and 



1 dp 



1 

< 

K T ds. 



cq 



cq 



(53) 



This means that the small thermal conduction, A and/or 
large relaxation time, r^j , is required in order to stabilize 
relativistic dissipative fluid. Recalling definition of the 
relaxation time for the heat fiow, = one of 

the stability conditions, ca. (l51ap . can be expressed as a 
constraint on the ^' appearing in eq. (j20b|) : 

l/e'<2(ecq + Pcq). 



Thus, conditions eg. ([55)1 for the stability of dissipative 
fluid can be expressed as the constraints for the inflnitcs- 
imal quantities 0(1) and 0(2) (See eq. ipOal) and (pOb)) ') 
in the following form: 



P 



cqj 



1 

^ < £cq 

1 dPcq 1 

^ < - < 1 



1 dp 



cq 



TdSe 



TdSe 



(54a) 
(54b) 



IV. 



SUMMARY AND CONCLUDING 
REMARKS 



We have proposed a new formulation of the relativistic 
dissipative hydrodynamical model in the Eckart frame by 
relaxing the standard matching conditions, 

Uf.uJT'"' = and u^SN^' = 0, 

and replacing by a more general form, 

Uf.uJT'"' = A and w^JiV^ = 6n. 

We assume that A = 0(1), Sn — 0(2) and that they are 
given by a simple function of dissipative quantities 11, 
W and tt'^^, cf. eq.^. 

Introducing the extended matching condition, we have 
also accordingly generalized the form of the entropy cur- 
rent for a non-equilibrium state in order that the ex- 
tended matching conditions might be suited (cf. eq. p^ ): 



S^q + W''/T- aSnui" + A/S" - x/^Hw^ 



The phenomenological parameter x introduced in the 
generalization of the entropy current form can be flxed 
by the extended thermodynamical stability condition (cf. 
eq-dlZl)) 

-fibSn -f A - xll = 0. 

As the result, the new the entropy current is formally 
the same as that obtained from Eckart theory (see 
eg. (lisp ). However, our matching conditions change 
the entropy production given by eg. (jl9p and therefore 
they also change the equation of motion of fluid (i.e., 
the constitutive equations 11, and tt'"^). This is 
because now x appears in the expression for the entropy 
production (cf. the 1st line of eg. dTOl) or, alternatively, 
in A and 6n cf. the 2nd line of eq. ([T9| ). It means then 
that, because of thermodynamic stability condition, the 
quantities of 0(1) and 0(2) introduced by the extended 



8 



conditions change both the entropy production and the 
equation of fluid motion. 

We check the stabihty of the fluid obtained by our 
model. Small perturbation imposed to a non-equilibrium 
state characterized by a constant bulk pressure Hq and 
= 0, TT^"^ = as the background reference fields. We 
found that the evolution of fluid is stable against small 
perturbations and the velocity of sound satisfies the usual 
limits, < Cs < 1, once eq. ()54a|) and eq. (j54b|) are sat- 
isfied. Those conditions give restrictions to the entropy 
production, eg. (1191) . via the constraints for k and ^' in 
the extended matching conditions eq. pOap and (j20bp . 

We conclude then, that when the matching condition 
Uf.uJT'"' = 0(1) and u^JiV^ = 0(2) are imposed prop- 
erly, the relativistic dissipative fluid described by the 
Eckart's first order theory evolves in stable and causal 
way. 



where x = 1 ~ k 



The third term is given by 



— TT^ 7 det Qa x 

+ [— + -hk') - Jlk' - -^/7eq] 



.AT 2 



J2fc2_|j^fc2 /icq 1^3 



-] 



+ [^(J3A:4 + _r^fc2)__£_J^fc2] p2 



4 , ^^2) 3f 
TRTtt {T) 

AT Js 



.AT 

■r, 

P AT 2 J3 ^ 4 Ji 



(A3) 



By adding eq. (|Al|) . eq. (|A2|) . and eq. (|A3|) . one obtains 
expression of eq. (|49|) . 



Appendix B: Sign of J2 
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J2 



de. 



oq 



ricq dT 



T dn, 



cq 
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consider the following thermodynamical relations 

rfPcq = ScqdT + UcqdH- (Bl) 



Appendix A: Determinant of the matrix Q 



The first term of eq.lgTl) times [j^^:^] is given by 



477 AT C 
(*fcr)det(Qi) X [-^ — ^] 

-^[^ + ^xk^]k^r^-^^^k^r, (Ai) 

^T^ TjjTw Tw OT-jY Tw TjJ 



where X = The second term is 



1-l-r^r 4?7 AT C 

— — — K ifc det Q2 X 

4?7/3 3t^ Tw Tu 

~ C r V , 1 ; 2-ri - C r ^P X 2a; 37^3 
— — K — [k X + x\k I — K — [K h — \k I 

Tj[ Tw Tjl Tw T-p; TwTi 

-f,^l^Xk^ + ^]k^r^-k^ — -k*T, (A2) 



Then, let's change of the thermodynamical variables from 
(ecq,"-oq) to (T, /i); 



dec 



C?77p 



dT 



^-dT- 



dn„ 



dT dfi 
With using a Jacobian defined by 



(B2a) 
(B2b) 

dT dfi dfi dT ' ^"^^^ 
one can inversely express relations eq. (|B2a|) and (jB2bp as 



J ^^cq ^^cq 



dn„ 



-dSn 



de, 



d^i ^^'^'^ 
-dsc 



dricq 



drio 
drir 



dfi ^ 

d, 



(B4a) 
(B4b) 



dT dT 
Substituting the above relations into eq. (|Bl|) . we have 



dPeq = 



Then we have 



J dfi 



J dT 



cq 



Scq dScq ^ 77cq dScq ^ 



J d^ 



J dT 



dPc 



cq 



Tn Tw 



Tn Tw TjT 



de, 



cq 



Scq dflcq 

J 9/i 



TT-cq dflcq 

~T~dT' 



(B5) 



(B6) 
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Eliminating /i from eq. (|B2ap and ea. (|B2bp . we have 
duca , de, 



and from this, one obtains 



dT 



9/i 



dn, 



cq 



eq 



1 de. 



JdT, 



cq 



(B7) 



(B8) 



Using eq. (jB6[) and (jB8|) . we have an expression 



:[Sc 



dUr, 



dT 



^cq d^cc\ 1 
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